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INTRODUCTION 

This  report  analyzes  Butterworth  filters  in  the  time  domain  by  using 
a  digital  method  for  finding  their  transient  performances.  Digital  filter 
compensation  is  effected  by  mviltiplication  in  the  Z- transform  domain  by 
(1+z™)/2.  For  the  sake  of  comparison,  posicast  compensation  (6)  of  a 
compensator  1  (1-P)+Pz"^  I  is  employed  in  each  corresponding  digital  filter 
compensation  case.  The  analytical  digital  process  used  in  this  report  is 
an  approximate  method  of  solving  differential  equations,  the  Naumov  program 
{k)t   which  is  an  improved  accuracy  trapezoidal  convolution  program  (3). 

Previous  work  on  Butterworth  filters,  pseudo-Butterworth  filters, 
posicast  control  and  Naumov  program  are  reviewed. 

Practical  examples  of  analyzing  and  improving  the  performances  of 
Butterworth  filters  from  second  order  to  fourth  order  with  unit  step  inputs 
are  presented. 

It  is  necessary  to  note  that  the  posicast  compensation  method  results 
in  larger  vmdershoots  and  time  delay  when  the  order  of  the  Butterworth 
filter  being  operated  becomes  higher.  These  defects  are  absent  in  digital 
filter  compensation.  Detailed  observations  are  presented  in  the  Discussion 
section, 

BUTTERWORTH  FUNCTIONS 
A  transfer  function  with  spectrum 

G^(juj)  =  -p=l=^  (1) 

yi+«/2n 


is  the  n^  order  Butterworth  filter  or  maximally  flat  low-pass  response  and 
it  is  an  approximation  to  the  ideal  low-pass  filter.  The  nature  of  the 
approximate  fimotion  is  seen  from  two  observations; 

(1)  From  the  binomial  series  expansion 

(l+X)-"=  XC^'^b  (-X)^     x2<  1  .  (2) 

k=0 

it  is  seen  that  the  Butterworth  response  when  UJ  is  near  zero  is 

(l  +  cO^r^/2  ^  i_(i/2)^2n  ^  (3/8)^^n  .  (5/16)0)^^  +  ....  (3) 
and  that  the  first  (2n-1 )  derivatives  are  zero  at  W=  0. 

(2)  [Sn^J'^^   ^*s  ^^^®<^  points  (  |g^^.J  »  |  G^  |  f or  all  n)  at 

u)=  0,  u)=   1  and  c^=  Oo   for  all  integers  n.  Corresponding  values  for 
Gn(Jt>0)   are  1,  0.70?,  and  0. 

The  Butterworth  transfer  functions  are  now  determined.  The  poles 
of  this  function  are  defined  by  the  equation, 

1  +  (-s2)^  =  0 
Then  the  pole  locations  are 

s,  =  e:5^2k-l)  n/Zn  „  even 

iv 

_  ^j(2k-hi-l);t/2n  k  =  1,2,3, t2n 

The  poles  so  defined  are  evenly  spaced  on  the  unit  circle  in  the 
s-plane  and  have  symmetry  with  respect  to  both  the  real  and  the  imaginary 
axes.  For  n  odd,  a  pair  of  poles  are  located  on  the  real  axis,  but  the 
poles  are  not  located  on  the  imaginary  axis  for  either  n  even  or  n  odd. 
These  properties  accrue  because  the  poles  are  separated  by  x/n  radians, 
and  axe  located  X/2n  radians  from  the  real  axis  for  n  even  and  on  the 
real  axis  for  n  odd.  ' 


si^  =  e^ 


or       ^k  ~  ® 


straightforward  calculation  of 

1 


Gn(s) 


aQ+ajS+a2s2+aos3+ +an^ 

_J w 


(s-sif)(s-Sk) 


yields  the  coefficients  aj^  of  the  denominator  polynomial  of  G^^Cs).  These 
Butterworth  polynomial  coefficients  are  given  in  Table  1  and  were  found 
in  reference  (?).  Here,  the  definition 

GnCs)  =  1/B„(s)  (M 

is  made* 

PSEUDO-BUTTERWORTH  FUNCTIONS 

The  pseudo-Butterworth  function  was  first  introduced  by  Marcel  J,  E, 
Golay  (2)  and  developed  by  K.  A.  Pullen,  Jr.  Pullen  (5)  extracted  the  basic 
polynomials  from  Golay 's  paper  and  worked  out  the  coefficient  relations. 

The  pseudo-Butterworth  polynomial  of  degree  n  is 

Bp(s)  =bo-H)^s-*2s2+ -*a-1^""^'^n^"  ^^^ 


=  Z:b.s^  (6) 


k=0 

This  polynomial  can  be  approximately  converted  into  a  true  Butter- 

l/n 
worth  type  by  replacing  s  with  s/(bjj)  '  .  The  approximate  Butterworth 

polynomial  is 

b'(s)  =  i:  a,  s^  .  •     (7) 

^     k=0  *^ 

where    a^^  =  hj{b^)^^''  (8) 


Table  1.  Coefficients  of  Butterwrorth  Polynomials 

Bn(s)  =  TT  (s-s^)(s-Si,) 

k=0 


I 


k 

0     1 


0  1.0 

1  1.0  1.0000 

2  1.0  \M^2  1.0000                                 .\  .     - 

3  1.0  2.0000  2.0000  1.0000                                .  ,• 
ii  1.0  2.6131  3.^1'^2  2.6131   1.0000 

5  1.0  3.2361  5.2361  5.2361   3.2361   1.0000 

6  1.0  3.8637  7M'^\  9.1^16   7.^^1   3.8637   1.0000       •' 

7  1.0  4.^^  10.0978  1^.5918  1^.5918  10.0978  >.^40   1.0000 

8  1.0  5.1258  13.1371  21.8462  25.688^  21.8462  13.1371   5.1258  1.0000 

The  coefficients  of  the  first  ten  orders  of  pseudo-Butterworth 
polynomial  and  their  corresponding  normalized  coefficients  of  the  first 
eight  orders  of  approximate  Butterworth  polynomial  are  presented  in  Table  2 
and  Table  3  respectively. 

The  basic  coefficient  recurrence  relation  of  an  n**^  order  pseudo- 
Butterworth  polynomial  is 

Cj  =  (n+j-2)  C^:]  +  Cj-''      0  ^  J  :^  n  (9) 

where  n  ^  2 

The  coefficients  of  zero  and  first  order  are 

cj  =  i,      c|=.i 

Note  that  the  superscripts  are  not  powers. 


Table  2.  Coefficients  of  Approximate  Butterworth  Polynomials 

B>)  =  g  s,."  .  a,  =  b./Cb^)"/" 

When  b,  's  are  the  coefficients  of  pseudo-Butterworth  polynomial 
k 


of  order 

n. 

1 

2 

? 

0  1.0 

1   1.0 

1 ,0000 

2  1.0 

1.^142 

1 .0000 

3  1.0 

2.0000 

2.0000 

1.0000 

^  1.0 

2.6596 

3.^87 

2.6392 

5  1.0 

3.3^5^ 

5.4633 

5.4388 

6  1.0 

if.  0^38 

8.0134 

9.9837 

1 .0000 

3.3684   1.0000 

7.8102   4.0000   1.0000 

7  1.0   4.7662  11.0079   16.0660   15.9415  10.7663   4.6125   1.0000 

8  1.0   5.4694  14.7004   25.2384   29.6392  24.8341  14.4298   5»4-718  1.0000 

The  additional  relations  are 

c"  =  1 
Lq  -  1 

C^=  (l/2)(n)^2^+1 

gJ  =  (l/8)(n+1)^^^+(l/2)(n+l)^2)+i 

C^  =  (l/48)(n+2)^^^+(l/8)(n+2)^^^-(l/2)(n-+2)^2^+3 

ej  =  (l/384)(n+3)^^^+(l/48)(n+3)^^^-(l/8)(n+3)^'*^+(3/2)(n+3)^^^-15 

C^  =  (l/3840)(n+4)^^°^+(l/384)(n-f4)^^^+(l/if8)(n+4)^^^+(3/8)(n-^)^'*'^ 

-(l5/2)(n+4)^2^+105 
etc. 
where  n^^^  f  n(n-1)(n-2) (n-k+1 ) 
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The  closed  form  for  C^  is 
^    2J  .  j  I      2^-''  .  ( j-1  )1   u=2         2n-^.(n-u)l       . 

and  the  pseudo-Butterworth  polynomial  can  be  written  as 

B  (s)  =  i  CV  (11) 


POSICAST  CONTROLS 

The  "posicast"  or  "positive-cast"  control  method  was  first  introduced 
by  Otto  J.  M,  Smith  (6),  TMs  method  can  be  applied  only  when  the  damping 
of  the  complex  poles  is  small.  The  secret  of  this  method  is  the  division 
of  the  input  function  into  two  times  and  two  amplitudes*  Ihe  second  input 
time  must  follow  the  first  input  time  l^y  exactly  one-half  cycle  of  the 
transient  if  the  superposition  of  the  results  is  to  be  zero  for  all  sub- 
sequent times.  Any  other  input-spacing  time  would  excite  two  oscillations 
whose  vectors  would  not  cancel,  because  they  would  always  have  some  phase- 
angle  difference.  Ihe  first  quantity  at  the  first  input  time  must  have  the 
value  that  can  bring  the  maximum  overshoot  to  reach  the  desired  final  out- 
put value  at  the  second  input  time.  The  second  quantity  must  have  the 
value  which  can  maintain  that  peak  value  produced  by  the  first  quantity. 

Figure  1  shows  a  natural  damped  transient  with  amplitude  A  at  one  peak 
and  amplitude  B  at  the  following  peak.  The   total  motion  is  A-fB.  Therefore 
a  step  A+B  should  be  broken  into  A  initially  and  B  later*  Ihe  control 
mechanism  should  have  the  transfer  function  as 


9 

p  ^  A  +  Be-^V2  ^^2) 

°   A  +B 

The  principle  of  vector  cancellation  with  posicast  control  can  also 

be  expressed  by  Jig.  1 .  In  Fig,  1 , 

(a)  Decrement  of  a  damped  oscillation.  ¥ector  A  which  starts  the 
oscillation  decreases  to  magnitude  B  after  time  T^/2»  (one-half  cycle 
later) . 

(b)  Required  excitation  function.  Vector  B  is  instantaneously  equal 
to  vector  A  and  180°  out  of  phase  with  it. 

(c)  Response  of  the  oscillatory  system  when  driven  by  the  excitation 
function  in  Fig.  1(b),  is 

A.u(t)  +  B»u(t  -  \i2) 
the  Laplace  transform  of  which  is  (A  +  Be"^  n/2)/s. 

r       ,      , 

THE  NAUMOV  PROGRAM 

Naumov  (^)  solved  differential  equations  in  the  time  dooain  employing 
trapezoidal  integration  on  a  multiple  integration  formula, 

d~y**  (t-^-)"  f(t)dr  (13) 

A  look-up  table  for  t^/nl  for  various  times  t  and  exponents  n  is  placed 
in  computer  storage.  A  method  is  then  given  for  employing  the  look-up 
table  in  conjunction  with  the  differential  equation  to  obtain  the  solution. 
Recently  Criswell  (1 )  devised  an  algorithm  that  easily  computes  exact 
Z(l/s^).  The  algorithm  readily  generates  tables  of  Z-transforms  of  t"/n! 
for  n  much  larger  than  those  in  presently  available  tables.  Such  a  large 
table  of  Z-transforms  can  be  used  for  checking  and  improving  approximate 


B 


(b) 


(c)    ■    •   .     ".    ^        • 
i"^g.l.  Principle  of  vector  cajicollation  with  Posicast  control 


10 

integration  formulas  generated  by  trapezoidal  convolution  (3).  Replacement 
of  approximate  Z- transforms  of  l/s^  by  their  exact  Z- transforms  yields  a 
trapezoidal  convolution  integrator  substitution  program  identical  to 
Naumov's  method  but  dispensing  with  digital  computer  memory  for  storing 
look-up  table  of  t^/nl  required  in  Naumov's  method,    •  ■ 

A  brief  derivation  of  Criswell's  algorithm  is  now  presented.  Since 
it  is  known  that 

Zd/s^)  -[xz/Cn-DJf^  Zd/s^-l)!  z  =  e-^=        (I'f) 

then  the  first  few  calculations  of  Z(l/s")  can  be  easily  carried  out. 
These  calculations  yield 

Z(l/s)  =  Tfz.sd  (15a)  ■ 

1   1-z 

Z(l/s2)  =Tz._l_^  (15b) 

1  (l-z)2 

Z(i/s3)=I^,11-t^  (150) 

2  (1  -  7.9 

Zd/s'*')  =  Ti,£!i5±l,         •  (15d) 

6   (1  -  z)^ 

_) 

Assume  that  the  previous  results  can  be  written  for  any  integer  n 
in  the  canonical  form 

(n-l)l   (1-z)'^ 

Where  Aj^(z)  is  a  polynomial  of  degree  n-2.  Differentiation  of  equation 
(l6)  shows  that  Z(l/s^  )  is  also  in  the  canonical  form.  Hence  by  mathematical 


11 

induction  the  canonical  form  holds  for  any  positive  n.  And,  the  differential- 
recurrence  relation  which  can  generate  the  desired  algorithm  is 

Aj^(z)  =  (1-z)  ^  (zAn_i)  +  (n-l)zAn_i  for  n  ^  2  (1?) 

Recall  that  the  starting  polynomial  is  A^(z)  =  z"' . 

Define  A(n,p)  to  be  the  coefficient  of  zP  in  A^^Cz).  The  differential- 
recurrence  relation  for  Aj^(z)  induces  the  coefficient  relation 

A(n,p)  =  (p  +  l)A(n-l,p)  +  (n-p-1  )A(n-1 ,p-l )  (18) 

The  coefficients  are  presented  in  Table  k  with  their  corresponding 
n's  and  p's. 

DIGITAL  FILTERS 

Conversion  of  a  linear  continuous  filter  function  into  digital  form 
can  be  accomplished  by  the  Naumov  program  in  the  Z-transform  domain.  15a.9 
linear  continuous  filter  functions  are  Butterworth  functions.  Since  the 
Naumov  program  gives  exact  Z(l/s^),  the  transient  response  of  a  digital 
filter  can  approach  the  exact  response  of  the  Butterworth  filter  by 
choosing  the  smaller  sampling  intejrval  sizes,  T.  The  numerical  calculations 
are  carried  out  by  a  digital  computer.  Conversion  of  the  Butterworth  filter 
into  a  digital  filter  with  unit  step  function  input  is  shown  for  the  second, 
third  and  fourth  order  Butterworth  filters,  „  - 

Second  Order  Digital  Butterworth  Filter 

If  the  second  order  Butterworth  filter  is 

GgCs)  = 1=:^ (19) 

^      1  +  >f2s  +  s2 
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then  the  output  function  with  \init  step  input  is 


- 

_  1 

1 

^ 

s     1 

+rzs 

+  s2 

)) 

=  y(o) 

=  0 

(20) 


2 
Rearranging  this  equation  and  dividing  tiy  s  gives 

0  +4I-  +  1    )y=  1  (21) 

s    s*^      s-* 
Taking  the  Z- transform  of  both  sides  yields 

Then,  by  employing  both  the  Naumov  program  and  the  approxuaate  Z-transform 

1  ** 
(3)  of  convolution,  (-^y).  gives 

ZC^y  )  =  '^(^+^)  Zy  (23a) 

2(1 -z)         ^  fs 

Z(J^)  =_l!^Zy  (23b) 


s' 


(1-z) 


,M)^T^z(1^)  ,  :.  ^^    ;.   •        (230) 

s3    2(1 -z)3  *■  ^-v  ^  ■   .; 

Substituting  equation  (23)  into  equation  (22)  and  then  multiplying  by 

2 
(1-z)  yields 


[z^d  -f)  +  z(t2  -  2)  +  (^  T  +  1)]  Zy 


_  T^z(1+z)  /pM 

2(i-z)  •  ::  ^^^^ 


»»^"B»J7r 
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The  corresponding  recurrence  relation  is 

y  =  2(2-1^)      +fe2       J^  (25) 


•  •  •  •  • 


X  =  jo,  1,  2,  2,  2,  2|  2,  •*•••••••••!  2, 

Let   A=2(2  -  T^) 
B=j2T  -  2 

D=|2T  -  2 
Then  y^=[(Ay^,l  +57^.2^/0]+  (C/D)x^  (26) 

When  n=0  yQ=0; 

n=l  yi=C/D; 

»^  V^^n-1  "■  ^n-2  -^  2C)/D. 

If  the  sampling  interval  T  is  chosen  as  0.1  second,  then  the  graphical 
and  numerical  solutions  of  y(t)  are  shown  in  fig.  2  and  Table  5  respectively. 
It  is  noted  that  the  graphical  solution  is  based  on  numerical  calculations. 

Third  Order  Digital  Butterworth  Filter 

The  third  order  Butterworth  filter  is 

B3(s)  = L-~ r-  (27) 

•^     1  +  2s  +  2s'=^  +  s^ 

Following  the  procedures  and  assumptions  of  the  second  order  case  yields 
Z  [(4^  +  ^  +  1 +  1)7J=  Z(Jjj.)  "    (28) 
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'able    5.   y(t)    =<j£."^-( 


lA 


s2+/23  +   1 


),y(0)    =0,    T   =0.1 


(t) 


0.0 

0.0000 

0.1 

0.0017 

0.2 

0.0180 

OoO 

0.0538 

0.4 

0.0685 

0.5 

0.0979 

0.6 

0.1345 

0.7 

0.1739 

0.8 

0.2160 

0.9 

0.2598 

1.0 

0.3048 

1.1 

0.3504 

1.2 

0.3959 

1.3 

0.4^'- 11 

1.4 

0.4856 

1.5 

0.5290 

1.6 

0.5710 

1.7 

0.6116 

1.8 

0.6513 

1.9 

0.6873 

2.0 

0.7222 

2.1 

0.7552 

2.2 

0,7861 

2.5 

0.8149 

2.4- 

0.8416 

2.5 

0.8663 

2.6 

0.8889 

y(t) 


2.7 

0.9096 

2.8 

0.9585 

2.9 

0.9455 

3.0 

0.9607 

3.1 

0.9743 

3.2 

0.9864 

3.3 

0.9970 

3.4 

1.0062 

3.5 

1 .0141 

3.6 

1.0209 

3.7 

1.0266 

3.8 

1.0313 

3.9 

lo0350 

4.0 

1.0380 

4.1 

1.0402 

4.2 

1.0417 

4.3 

1.0426 

4.4 

lo0431 

4o5 

1,0430 

4.6 

1.0426 

4.7 

1,0418 

4o8 

1.0408 

4.9 

1.0394 

5.0 

1,0379 

5,1 

1.0563 

5,2 

1,0345 

5o3 

1,0326 

5.5 

5.6 

5.7 

5.8 

5.9 

6.0 

6.1 

6.2 

6.3 

6.4 

6.5 

6.6 

6.7 

6.8 

6.9 

7.0 

7.1 

7.2 

7.3 

7.4 

7.5 

7.6 

7.7 

7.8 

7.9 


y(t) 


1.0307 

1,0288 

1 .0268 

1.0248 

1.0229 

1.0210 

1.0192 

1.0174 

1.0157 

1.0141 

1.0126 

1.0111 

1.0097 

1.0085 

1.0073 

1.0062 

1.0052 

1.0042 

1.0034 

1.0026 

1.0019 

1.0013 

1.0008 

1.0003 

0.9999 

0.9995 
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The  Naumov  program  gives 

K  _  T^z  z  Vfz+1  (29d) 


Then 


Zy  [(T-1 )z3+(3-T-2t240.5T^)z2+(0.5T^+2t2-T-3)z+(T+1  )  j 

=  £sl£±ih5+li  (30) 

6         (1-z) 


Let 


A  =  (3+T-2T^-0.5T^) 

B  =  (t+2t2-o.5t3-3)  '. 

C  =  1-T  ■  , 

D  =  T^ 
R  =  1+T 
Then  the  recurrence  relation  is 

yn=[(^n.l^^n-2-^n-3^/^J^^^^^      '/.     ^^'^ 


When 


X  =  Jo,  If  5t  ^>  ^t  6*  •••••••  o»  ••••  > 

n  =  0      70  ""  °' 

n  =  2     yz'I'l  *§' 
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n^3      yn=^^n-l'"^n-2'^n-3°^/^- 
Results  are  presented  in  Jig.  2  and  Table  6, 

Fourth  Order  Digital  Butterworth  Filter 
The  fourth  order  Butterworth  filter  is 

^^        "  s^+2.6s3+3.^s2+2.6s+1  (32) 

Calculation  yields  ■,  ,  .     - 

zfcl  +2^+M.^h-1)jf]=Z(1)  (33) 

U  s^  s-'   s^    s 

The  Navmov  program  gives  : .  . 

7(2. 6y)  =    1.3T(1+z)   7y  (3^a) 

s2  (1-z) 

Z(My)  =   Iiii£S-Zy  (3^) 

s3   "^  (1-z)2 

Z(2iiy)   =    lLl£sll+2Lzy  (3^) 

83  (1-2)3 

z(Jpy)  =  l!:si£fi:5+il-zi^  r  (3iki) 

s^  6(1-z)^ 

„(1  )  -  A(z3+llz2+llz-H)        '  (3i^) 

s5  2Ml-z)5 

Let  W  =  ^ 

o 


E  =  2- 


i9 


Table 

6.    y(t)   =X 

'■(  .^sK 

2s^+2s+l)     ' 

,  y(o)=o. 

T=0.1    sec. 

t 

y(t) 

t 

y(t) 

t 

y(t) 

0.0 

0.0000 

3.2 

0.8750 

6.4 

1.0300 

0.1 

0.0002 

3.3 

0.9012 

6.5 

1.0256 

0.2 

0.0012 

3.4 

0.9254 

6.6 

1.0214 

0.3 

0.0038 

3.5 

0.9477 

6.7 

1.0174 

0.4 

0.0086 

3.6 

0.9681 

6.8 

1.0137 

0.5 

0.0160 

3.7 

0.9866 

6.9 

1.0101 

0.6 

0.0263 

3.8 

1.0033 

7.0 

1.0067 

0.7 

0.0397 

3.9 

1.0181 

7.1 

1.0036 

0.8 

0.0562 

4.0 

1.0311 

7.2 

1.0008 

0.9 

0.0758 

4.1 

1.04  24 

7.3 

0.9982 

1.0 

0.0985 

4.2 

1.0521 

7.4 

0.9959 

1.1 

0.1242 

4.3 

1.0602 

7.5 

0.9938 

1.2 

0.1525 

4.4 

1.0668 

7.6 

0.9920 

1.3 

O.I834 

4.5 

1.0721 

7.7 

0.9904 

l.h 

0.2165 

4.6 

1.0760 

7.8 

0.9891 

1.5 

0.2515 

4.7 

1.0788 

7.9 

0.9890 

1.6 

0.2882 

4.8 

I.O8O5 

8.0 

0.9871 

1.7 

0.3263 

4.9 

1.0811 

8.1 

0.9864 

1.8 

0.3654 

5.0 

I.O8O9 

8.2 

0.9859 

1.9 

0.4053 

5.1 

1.0798 

8.3 

0.9856 

2.0 

0.4465 

5.2 

I.O78I 

8.4 

0.9854 

2.1 

0.4861 

5.3 

1.0757 

8.5 

0.9854 

2.2 

0.5264 

5.4 

1.0727 

8.6 

0.9856 

2.3 

0.5663 

5.5 

1.0693 

8.7 

0.9858 

2.k 

0.6056 

5.6 

1.0656 

8.8 

0.9862 

2.5 

0.6441 

5.7 

1.0615 

8.9 

0.9867 

2.6 

0.6815 

5.8 

1.0572 

9.0 

0.9872 

2.7 

0.7176 

5.9 

1.0528 

9.1 

0.9878 

2.8 

0.7524 

6.0 

1.0482 

9.2 

0.9885 

2.9 

0.7856 

6.1 

1.0436 

9.3 

0.9892 

3.0 

0.8172 

6.2 

1.0390 

9.4 

0.9900 

3.1 

0.8470 

6.3 

1.0344 

9.5 

0.9907 
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Substituting  equation  (3^)  into  equation  (33)  yields 

•      z7r(1+^.3T)+z(-i>-2.6T+3.^T2+l.3T^-*W)  +  z2(6-6.8t2+W)  + 

z3(-i^+2.6T+3.^T2-U3T^-*W)  +z^(l-1.3T)j 
E5(z3-;-11z^+11z+1)  (35) 

(l-a)  .;■■•        .     ::■•/:    . 

Let  A  =  ^+2.6T-3.^T^-1.3T^-W 

B  =  6.8T^-W-6  "  ■-■■   ■    \-   '.^  ' 

C  =  ^-2.6T-3.^T^+1  .3T^-W  '.  .  '■  |\ 

D  =  1.3T-1 
F  =  1+1.31 


Then 


^n 


F                                ■  F         . 
x^^  =  jo,   1,   12,  23,  2ii',  24,   24,    24 1- 

■VKien  n  =  0  ^n  "  *^'   ' 

•      n  =  2  72  =  (Ay^+12E)/F; 

n  =  3  y^  =  (A72+Bjr^+23E)/F: 

n  ^  4  y^  =  (Ay^.i-^By^.2^^.3-*-^^^+24E)/F. 

Results  are  sbovm  both  in  Table  7  a-nd  Fig.  2.  "■  /  " ;- 
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Table  ?.     y(t)  =  ^"^  ( 


111 


s^+2 . 6  s3+3 . /j.s2+2 . 6  s+1 


-)t  y(0)=0,  T=0.1  seo. 


y(t) 


0.0 
0.1 
0.2 

0.3 

O.ij- 

0.5 
0.6 

0.7 

0.8 

0.9 
1.0 
1.1 
1.2 
1.3 
1.4 

1.5 
1.6 

1.7 
1.8 

1.9 
2.0 
2.1 
2,2 
2.3 

2.5 
2.6 

2.7 
2.8 

2.9 
3.0 
3.1 
3.2 

3.3 
3.4 

3.5 
3.6 

3.7 

3.8 
3.9 
4.0 
4.1 
4.2 

4.3 

4.4 


0.0000 
0.0000 
0.0000 
0.0003 
0.0009 
0.0020 

0.0039 

0.0068 

0.0111 

0.0168 

0.0242 

0.0334 

0.0447 

0.0581 

0.0738 

0.0917 
0.1118 
0.1342 
0.1588 

0.1855 
0.2142 
0.2447 
0.2770 
0.3107 
0.3457 
0.3818 
0.4189 
0.4565 
0.4946 
0.5330 
0.5713 
0.6093 
0.6469 
0.6839 
0.7200 

0.7551 
0.7891 
0.8217 
0.8528 
0.8824 
0.9103 
0.9365 
O.9608 

0.9833 
1 .0039 


4.5 
4.6 

4.7 
4.8 
4.9 
5.0 
5.1 
5.2 
5.3 
5.4 

5.6 

5.8 
5.9 
6.0 
6.1 
6.2 
6.3 
6.4 

6.5 
6.6 

6.7 
6.8 

6.9 
7.0 
7.1 
7.2 
7.3 
7.4 

7.5 
7.6 
7.7 
7.8 

7.9 
8.0 
8.1 
8.2 

8.3 
8.4 

8.5 
8.6 

8.7 

8.8 

8.9 


y(t) 

t 

y(t) 

1 .0226 

9.0 

0.9714 

1 .0394 

9.1 

0.9717 

1 .0543 

9.2 

0.9723 

1 .0673 

9.3 

0.9732 

1 .0785 

9.4 

0.9742 

1 .0880 

9.5 

0.9754 

1 .0957 

9.6 

0.9768 

1.1019 

9.7 

0.9784 

1.1064  • 

9.8  . 

.  '   0.9800 

1.1096 

9.9 

O.98I7 

1.1113     1 

0.0 

0.9835 

1.1118     1 

0.1 

.  0.9854 

1.1111     1 

0.2 

0.9872 

1.1093     1 

0.3 

0.9891 

1.1066     1 

0.4 

0.9909  . 

1.1030     1 

0.5 

0.9927 

1 .0987     1 

0.6 

0.9954 

1 .0937     1 

0.7 

0.9962 

1 .0882     1 

0.8 

0.9978 

1 .0822 

0.9 

0.9994 

.  1.0759 

1.0 

1 .0008 

1 .0693 

1.1 

1 .0022 

1 .0626 

11.2 

1.0034 

1.0557 

11.3 

1 .0045 

1 .0488 

1.4 

1 .0056 

1 .0420 

11.5 

1.0065 

1 .0353 

11.6 

1 .0073 

.  1 .0287 

11.7 

1 .0079 

1 .0224 

11.8  ■ 

1 .0085 

1.0163 

11.9 

1 .0090 

■  1.0106 

12.0 

1 .0093 

1 .0052 

12.1 

■  1 .0095 

1 .0001 

12.2 

1 .0097 

0.9955 

12.3 

1 .0098 

0.9912 

12.4 

1 .0097 

0.9874 

12.5 

1 .0096 

0.9840 

12.6 

1.0094 

0.9810 

12.7 

1.0092 

0.9785 

12.8 

1 .0089 

0.9763 

12.9 

1.0085 

0.9746 

13.0 

1.0081 

0.9732 

13.1 

1 .0077 

0.9723 

13.2 

1 .0072 

0.9717 

13.3 

1 .0068 

0.9714 

13.4 

1 .0063 
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DIGITAL  FILTER  COMPENSATION 

Posicast  control  compensation  can  only  be  applied  when  the  damping 
is  small  and  introduces  large  time  delay  in  the  transient  response;  as 
the  order  of  the  system  becomes  higher,  posicast  compensation  always 
results  in  large  undershoots  in  the  transient  response.  For  the  sake 
of  reducing  these  defects,  "digital  filter  compensation"  is  employed. 

The  basic  idea  of  digital  filter  compensation  is  conversion  the 
unit  step  function  input  into  the  double  half-step  function  displayed 
in  Fig.  3(b),  The  second  half-step  input  follows  the  first  half-step 
function  after  a  time  mT  seconds,  where  mT  is  the  time  duration  of  the 
first  overshoot  of  the  Butterworth  response  to  a  unit  step  function.  This 
is  clarified  in  Rig,  3(a). 

The  symbol  T  is  the  duration  between  sampling  times,  and  m  is  a 
constant  to  be  determined  for  each  Butterworth  filter.  Therefore  a 
pre-compensating  operater 

1  -  e"°^^^  =  1  -  z°^  (37) 

2  2 

can  be  constructed.  The  block  diagram  of  the  system  with  this  compensator 

is  shown  in  Fig.  7(a), 

Where  X,  Y  and  G(s)  are  the  input,  output  and  transfer  (Butterworth 

filter)  functions  respectively.  This  completes  the  first  step  of  compensation. 

This  step  reduces  most  of  the  overshoot  of  the  original  transient  response 

but  as  happens  in  posicast  control  compensation,  there  now  exists  larger 

time  delay  which  is  much  greater  than  that  in  the  original  system.  Hence, 

the  next  step  is  to  induce  time  delay. 
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M 


1.0 


0.5 


t  a  mT 


(b) 


Pi^.  3.  Principle  of  double  atep  input  function-  of  digital 
■  filter  compensation 


2^ 

Assume  that  the  delay  time  in  the  original  system  is  T^^  and  after 

first  step  compensation  is  T^.  Since  T^  is  greater  than  1^^,   the  ratio 

T,  /T  ,  is  some  value  less  than  1 .  The  procedure  is  to  replace  s  by 
do  dl 

(Tj  /Tj,)s  in  the  compensated  Butterworth  polynomial.  This  replacement 
do  dl 

of  s  will  renormalize  time  delay  since  T^^/T^  is  less  than  1 ,  this  meets 
the  second  requirement;  so,  the  goal  of  compensation  Is  accomplished. 

Transient  behaviors  of  several  orders  of  Butterworth  filter  by 
applying  digital  filter  compensation  will  be  investigated  in  the  next 
section. 

Compensation  of  Second  Order  Digital  Butterworth  Filter 

A,  First  step  of  compensation  '■,.'. 

The  time  duration  mT  of  the  first  overshoot  of  G2C2)  with  unit  step 
input  should  be  known;  this  can  be  found  from  Table  5  where 
T  =  0.1  sec, 
mT  =  7,8-3.^  =  ^A   sec. 
Therefore  m  =  ^. 

Then  the  system  response  with  unit  step  input  is  shown  in  Pig.  7W* 
Thus  .    1  ^  ^ 

y  =  0-)  (-i-f^)(      1        )    .  -  (38) 

^     2     s2  +  V2s  +  1 

1   1  +  z^ 
Rearranging  equation  (38)  by  leaving  (-)( — r )  at  one  side  yields 

(s2+/2s+l)y=(l)Li.^  .     (39) 

2 

Taking  the  Z-transform  yields  after  division  by  s  at  both  sides 

After  applying  the  Naumov  program  this  equation  yields 
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Zy  [(1  -  £t)z2  +  (t2-2)z  +  (1  +^T)] 


i^(l-z) 


(^1) 


Let  A  =  2(2-T^) 

B  =4^1-2 

C  =  T^ 

D  =-f2T+2  '  • 

Then  the  recurrence  relation  is  ^ 

y  ,  (^n-1^Syn-2)  ^(^)  '(,,2) 

X  =  0,  1 ,  2|  ••«•  2,  3»  ^»  ^»  •••••  ^t    ••••• 

This  equation  yields 

y  =0  n  =  0 

y  =^  n  =  1 

'2D 

^n  =  (^n-l''^n-2'^^/°  2  ^  n  <44 

y45  =  (Ay/^4+B7^3+1.5C)/D  n  =  ^5 

^n  "  ^^n-l'^^n-2^^^^/°  n  ^46  ' 

The  numerical  and  graphical  results  of  y(t)  are  presented  in  Table  8 
and  Fig.  ^  respectively, 

B,  The  second  step  of  compensation 

From  Table  5  and  Table  8,  the  delay  times  which  are  found  before  and 
after  the  first  step  compensation  are 

■^do  ~  "*  '^^   ^®°' 
"^dl  ^  ^'^^  ^^^' 
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Table    8.    y(*)  =  f''''( 


1  +  z 


kh 


lli 


■)»    y(0)=0»    T=0.1   sec. 


2s    +    1 


y(t) 


0.0 

0.0000 

0.1 

0.0023 

0.2 

0.0090 

0.3 

0.0194 

0.4 

0.0329 

0.5 

0.0490 

0.6 

0.0671 

0.7 

0.0896 

0.8 

0.1080 

0.9 

0.1299 

1.0 

0.1524 

1.1 

0.1752 

1.2 

O.I98O 

1.3 

0.2206 

1.4 

0.2428 

1.5 

0.2645 

1.6 

0.2855 

1.7 

0.3056 

1.8 

0.3252 

1.9 

0.3436 

2.0 

0.3611 

2.1 

0.3776 

2.2 

0.3930 

2.3 

0.4074 

2.k 

0.4208 

2.5 

0.4331 

2.6 

0.4445 

2.7 

0.4548 

2.8 

0.4642 

2,9 

0.4727 

3.0 

0.4804 

3.1 

0.4872 

3.2 

0.4932. 

3.3 

0.4985 

3.4 

0.5031 

3.5 

0.5071 

3.6 

0.5105 

3.7 

0.5133 

3.8 

0.5156 

3.9 

0.5175 

4.0 

0.5190 

4.1 

0.5201 

4.2 

4.3 
4.4 

4.5 
4.6 

4.7 
4.8 

4.9 
5.0 

5.1 
5.2 

5.3 

5.4 
5.5 
5.6 
5.7 
5.8 

5.9 
6.0 
6.1 
6.2 

6.3 

6.4 

6.5 
6,6 

6.1 
6.8 

6.9 
7.0 

7.1 
7.2 

7.3 

7.4 

7.5 
7.6 

7.7 
7.8 

7.9 
8.0 
8.1 
8.2 
8.3 


y(t) 


0.5209 
0.5213 
0.5215 
0.5239 
0.5303 
0.5403 
0.5533 
0.5687 
0.5861 

0.6051 
0.6252 
0.6462 
0.6678 
0.6896 
0.7114 
0.7330 
0.7543 
0.7750 
0.7951 
0^8145 
0.8330 
0.8507 
0.8674 
0.8831 
0.8979 
0.9117 
0.9244 
0.9362 
0.9470 
0.9569 
0.9659 
0.974b 
0.9813 
0.9878 
0.9936 
0.9986 
1.0030 
1.0068 
1.0101 
1.0128 
1.0150 
1.0168 


8.5 

8.6 

8.7 
8.8 

8.9 
9.0 

9.1 

9.2 

9.3 
9.4 
9.5 
9.6 
9.7 
9.8 

9-9 
10.0 
10.1 
10.2 

10.3 
10.4 

10.5 
10.6 

10,7 
10.8 

10.9 
11.0 
11.1 

11.2 

11-3 
11.4 

11.5 
11.6 

11.7 
11.8 

11.9 
12.0 
12.1 
12.2 

12.3 
12.4 

12.5 


(t) 


1.0192 
1.0200 
1.0204 
1.0206 
1.0206 
1.0204 
1.0200 
1.0195 
1.0189 
1.0181 
1.0173 
1.0165 
1.0156 
1.0147 
1.0137 
1.0128 

1.0119 
1.0109 
1.0100 
1.0092 
1,0083 
1.0075 
1.0067 
1.0060 

1.0053 

1,0046 

1.0040 

1.0035 

1.0029 

1.0025 

1.0020 

1.0016 

1.0013 

1.0009 

1.0006 

1.0004 

1.0001 

0.9999 

0.9998 

0.9996 

0.9995 
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Then 


_do  =  Iti]:!  =  0.-4-27 
^dl   3.35 


Replacing  s  by  0.^27s  in  the  transfer  function  (  ^ 
of  equation  (38)  yields 


■  y  ~  ^s^^       2    '^S 40.605+0.01 81  s2 


1       \ 

1  +//rs  +  s^ 


(^3) 


Rearranging  and  tald.ng  the  Z-transform  of  both  sides  yields 

Z  [(-L  +  0^  +  0.0l81)y]  =  z[^(l-|-2li)j    C^) 

Application  of  the  Naumov  program  yields 

Zy  [(0.0181  +  0.303T)  +  (t2-0.0362)z  +  (0.0181-0.303T)z2j 

=  T^2(1+a)(T+z^^)  (^5) 

4(1-z) 

Let  A  =  0.0362-t2 

B  =  0.303T-0.0181 

C  =  T^ 

D  =  0.303T  +  0.0131 
Then,  the  recurrence  relation  is 

_(^y^-1-^By^-1)  -.(0.25CW 

Xj^  =-fo,  1,  2,  2 2,  3,  ^,  ^. 

This  ..equation  yields 

,  yo  =  ° 

y^    =  0.25C/D 

yn  =  ^^n-l-^By^-2^.5C)/D 

^20  =  (Ayi9+Byi8+0'75C)/D 


(^6) 


^1 

>  •  •  .  I    "r    r 

n  =  0 
n  =  1 
2  :^n  ^19 

n  =  20 
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y^=(Ay^.l+By^.2^)/D  n  ^  21 

The  graphical  and  numerical  results  are  shown  in  Fig.  ^  and  Table  9 
respectively.  In  Jig.  4,  we  can  see  that  the  transient  response  after  the 
second  step  of  compensation  is  better  than  the  original  Butterworth  filter 
response.  The  maximum  overshoot  is  reduced  from  0.0^31  to  none  but  the 
delay  time  increases  to  0.52  seconds. 

Compensation  of  Third  Order  Digital  Butter.Jorth  Filter 

A,  RLrst  step  of  compensation 
Table  6  gives  m  =  35 
Thus 

The  block  diagram  of  this  system  is  shown  in  Fig.  7(o). 
Calculation  yields 

Zy  [z^CT-D  +  z2(3-T-2t2+0.5T^)  +  z(0.5T^-^2T2-T-3)  +  (T+1)] 

^:-         =  T^zCz^-f^z+Dd+z^^)  "  (48) 

I2(l-z) 

Then  the  recurrence  relation  is 

-  (Ay^.i4-By„.2->-Cyn_3)     ,    p  _„  V''       •  ^^^  . 

^  a  12R^  . 

x^  ~'|^t  1>  5i  6f  ^t    .....>  6,  7t   11t  12,  12,   ......  12,   .....  y 

A  =  3+T-2t2-0.5T^ 
B  =  T+2t2-0.5t5-3 
C  =  1-T 
D  =  t3 
R  =  1+T 


Where 
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9.  y(t)=X" 

'^   2 

1/s 

-  '\ 

Table 

/   ,    ,2 

2  '* 

t+  2(0.427)s+ 

(0.427) 

s 

y(o)= 

0,  T=0.1  sec. 

t 

y(t) 

t 

y(t) 

t 

y(t) 

0.0 

0,0000 

2.7 

0.8634 

5.4 

0.9987 

0.1 

0.0025 

2.8 

0.8851 

5.5 

0.9989 

0.2 

0.1047 

2.9 

0.9034 

5.6 

0.9991 

0.3 

0.1606  ■ 

3.0 

0.9188 

5.7 

0,9992 

o.i; 

0.2166 

3.1 

0.9317 

5.8 

0.9994 

0.5 

0.2610 

3.2 

0.9425 

5.9 

0,9995 

0.6 

0.2992 

3.3 

0.9517 

6.0 

0.9995 

0.7 

0.3311 

3.4 

0.9594 

6.1 

0.9996 

0.8 

0.3579 

3.5 

0.9658 

6.2 

0.9997 

0.9 

0.3805 

3.6 

0.9712 

6.3 

0.9997 

1.0 

0.3995 

3.7 

0.9758 

6.4 

0.9998 

1.1 

0.4155 

3.8 

0.9797 

6.5 

0.9998 

1.2 

0.4289 

3.9 

0.9829 

6.6 

0.9999 

1.3 

0.4402 

4.0 

0.9856   :"  = 

6.7 

0.9999 

1.4 

0.4497 

4.1 

0.9B79 

6.8 

0.9999 

1.5 

0.4577 

4.2 

0.9898 

6.9 

0.9999 

1.6 

0.4644 

4.3 

0.9914 

7.0 

0.9999 

1.7 

0.4701 

4.4 

0,9928 

7.1 

0.9999 

1.8 

0.4748 

4.5 

0.9939 

7.2 

0,9999 

1.9 

0.4788 

4.6 

0.9949 

7.3 

1,0000 

2.0 

0.5339 

4.7 

0.9957 

7.^ 

1,0000 

2.1 

0.6163 

4.8 

0.9964 

7.5 

1.0000 

2.2 

0.6748 

4.9 

0.9970 

7.6 

1.0000 

2.3 

0.7272 

5.0 

0.9975 

7.7 

1.0000 

2.k 

0.7704 

5.1 

0.9979 

7.8 

1.0000 

2.5 

O.8069 

5.-2 

0.9982 

7.9 

1.0000 

2.6 

0.8376 

5.3 

0.9985 

8.0 

1,0000 

31 


Equation  (^9)  yields 


y     =  0  n  =  0 


Then 


T^  =  3.75     sec. 


-do  =  2JiZ  =  0.573 
Tdi      3.75 


n  =  1 
n  =  2 


0 
y^  =  D/12R 

y^  =  (Ayi/R)  +  (5D/12R) 

y^  =  (Ay^_^  +By„^2+Cy„.3-«5-5D)/R  3  <  n  4  35 

^36  =    CAy33+By3i^-K:y33)/R    +  7R/12R  n  =  36 

y37  =    (Ay36+^735+<^3i^)/R    +  11D/12R  n  =  37 

yn=  Uyn-l+Byn.2^n-3"*"°^/^  ""^^^ 
Results  are  presented  in  Table  10  and  Fig.  5. 

B.  Second  step  of  compensation 

From  Table  6  and  Table  10,  we  get 

■^do  ^  ^•'•^'^  ^®°* 


Replacing  s  by  0.573s  in  the  transfer  function  of  equation  (^7)  yields 


y  =  (s)(— 2    ^^0.188s3-K).56is2+1.146s+l  ' 


(50) 


Applying  the  same  process  as  in  second  order  case  yields 

ZyRo. 188+0. 327T)+(1.1^6T^-0.327T-0. 56^+0. 5T^)z 
+(0.106-0.327T-1.1^6t2+0.5t3)z2+(o.327T-0.188)z3J 

12(1-z) 
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Table  10.  y(t)=/~  (~^ 


.35 


JJJl 


s-^+2s^+2s+l 


),  y(0)=0,  T=0.1  sec. 


y(t) 


0.0 

0.0000 

0.1 

0.0000 

0^2 

0.0006 

0.3 

0.0019 

o.k 

0.0043 

0.5 

0.0080 

0.6 

0.0132 

0.7 

0.0199 

0.8 

0.0281 

0.9 

0.0379 

1.0 

0.0432 

1.1 

0.0621 

1.2 

0.0763 

1.3 

0.0917 

1.4 

0.1082 

1.5 

0.1258 

1.6 

0.1441 

1.7 

0.1632 

1.8 

0.1827 

1.9 

0.2027 

2.0 

0.2228 

2.1 

0.2430 

2.2 

0.2632 

2.3 

0.2832 

2,h 

0.3028 

2.5 

0.3220 

2.6 

0.3407 

2.7 

0.3588 

2.8 

0.3762 

2.9 

0.3928 

3.0 

0.4086 

3.1 

0.4235 

3.2 

0.4375 

3.3 

0.4506 

3.^+ 

0.4627 

3.5 

0.4739 

3.6 

0.4841 

3.7 

0.4939 

3.8 

0.5035 

3.9 

0.5134 

4.0 

0.5236 

4.1 

0.5344 

4.2 

0.5459 

4.3 

0.5582 

4.4 

0.5713 

4.5 

4.6 

4.7 
4.8 

4.9 
5.0 

5.1 
5.2 

5.3 
5.4 
5.5 
5.6 
5.7 
5.8 

5.9 
6.0 
6.1 
6.2 

6.3 

6.4 

6.5 
6,6 

6.7 
6.8 

6.9 
7.0 

7.1 
7.2 

7.3 
7.4 

7.5 
7.6 

7.7 
7.8 

'7.9 
8.0 
8.1 
8.2 

8.3 
8.4 

8.5 
8.6 

8.7 
8.8 

8.9 


y(t) 


y(t) 


0.5853 
0.6001 
0.6157 
0.6319 
0.6488 
0.6662 
0.6840 
0.7022 
0.7206 
0.7390 
0.7575 
0.7758 

0.7939 

0.8118 

0.8292 

0.8461 

0.8625 

0.8783 

0.8934 

0.9078 

0.9214 

0.9342 

0.9462 

0.9574 

0.9677 

0.9772 

0.9859 

0.9937 

1.0007 

1.0070 

1.0125 

1.0172 

1.0212 

1.0246 

1.0274 

1.0296 

1.0312 

1.0323 

1.0330 

1.0333 

1.0331 

1.0327 

1.0319 

1.0309 

1.0297 


9.0 
9.1 
9.2 
9.3 
9-4 

9.5 
9.6 

9.7 
9.8 

9.9 
19.0 
10.1 
10.2 

10.3 
10.4 

10:5 
10.6 

10.7 
10.8 

10.9 
11.0 
11.1 
11.2 

11.3 
11.4 

11.5 
11.6 

11.7 
11.8 

11.9 
12.0 
12.1 
12.2 

12.3 
12.4 

12.5 
12.6 

12.7 
12:8 

12.9 
13.0 

13.1 
13.2 
13.3 
13.4 


1.0282 

1.0267 
1.0250 
1.0232 
1.0213 
1.0194 

1.0175 
1.0156 

1.0137 
1.0118 
1.0101 
1.0083 
1.0067 
1.0051 
1.0037 
1.0023 
1.0010 

0.9999 
0.9988 

0.9979 
0.9970 

0.9963 
0.9957 
0.9952 

0.9947 
0.9944 

0.9941 
0.9939 
0.9938 
0.9938 
0,9938 
0.9939 
0,9940 
0.9942 
0,9944 
0.9946 

0.9949 
0.9952 
0.9955 
0.9958 
0.9962 

0.9965 
0.9968 
0.9972 
0.9975 
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Then  the  recurrence  relation  is 

rAyn..1-^By^-2^n-3)     .  DX.  (52) 

"  R  12R 

x^  =fo,1,5,6,6 .6,7,11,12,12 12 I 

Where 

A  »  0.56'4- +  o.327T-i.i^t^-o.5t2 

B  =  O.327T  +  1.1^6T^-0.5T^-0.56^ 

C  =  0.1 88-0. 327T 
;  .         D  =  "P 

R  =  0.188  +  O.327T  .... 

Equation  (52)  yields 

yo  =  °  '■; 

71  =  D/12R 

72  =  (Ay^/ft)  +  (5D/-12R) 
yn  =  ^^yn-l^^n-2'^n-3'^*^°^/^  3  ^  n  ^20 

721  =     (A720+S7i9+C7i8)/^     +(7D/12R) 

722  =     (Ay2i+B72o+C7i^)/a     +(11D/12R) 

yn  =  (^n-l'"^n-2'*^n-3"'°^/^  ""  ^^^ 

The  graphical  and  numerical  results  are  presented  in  Fig.  5  and 
Table  11  respectively.  The  maximum  overshoot  reduces  from  0.0811  to 
0.0328  and  the  time  dela7  reduces  from  2.157  seconds  to  2.15  seconds. 

Compensation  of  Fourth  Order  Digital  Butterworth  Filter 

A.  First  step  of  compensation 
Table  7  gives  m  =  33. 


n  =  0 
n  =  1 
n  =  2 


n  =  21 
n  =  22 
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Table  11.  y(t)=o£'"-'-(  — 


20 


(t) 


^    0.l88s^+0.65^s^+l.lit6s+l 
y(0)=0,  T=0.1  sec. 


). 


y(t) 


0.0 

0.0000 

4.2 

1.0033 

0.1 

0.0004 

4.3 

1.0131 

0.2 

0.0029 

4.4 

1.0207 

0,3 

0.0090 

4.5 

1.0262 

0.4 

0.0196 

4.6 

1.0299 

0.5 

0.0350 

4.7 

1.0320 

0.6 

0.0552 

4.8 

1.0328 

0.7 

0.0796 

4.9 

1.0324 

0.8 

0.1078 

5.0 

1.0310 

0.9 

0.1390 

5.1 

1.0289 

1.0 

0.1725 

5.2 

1.0263 

1.1 

0.2072 

5.3 

1.0233 

1,2 

0.2426 

5.4 

1.0201 

1.3 

0.2778 

5.5 

1.0168 

1.4 

0.3120 

5^6 

1.0135 

1-5 

0.3449 

5.7 

1.0104 

1.6 

0.3757 

5.8 

1.0074 

1.7 

0.4043 

5.9 

1.0047 

1.8 

0.4301 

6.0 

1.0023 

1.9 

0.4532 

6.1 

1.0002 

2.0 

0.4733 

6.2 

0.9984 

2.1 

0.4909 

6.3 

0.9969 

2.2 

0.5078 

6.4 

0.9958 

2.3 

0.5255 

6.5 

0.9949 

2.4 

0.5451 

6.6 

0.994^ 

2.5 

0.5671 

6.7 

0.9940 

2.6 

0.5917 

6.8 

0.9939 

2.7 

O.6187 

6.9 

0.9940 

2.8 

0.6478 

7.0 

0.9942 

2.9 

0.6785 

7.1 

0.9946 

3.0 

0.7103 

7.2 

0.9950 

3.1 

0.7427 

7.3 

0.9956 

3.2 

0.7749 

7.4 

0.9961 

3.3 

0o8064 

7.5 

0.9967 

3.4 

0.8369 

7.6 

0.9973 

3.5 

0.8657 

7.7 

0.9978 

3.6 

0.8926 

7.8 

0.9983 

3.7 

0.9173 

7.9 

0.9988 

3.8 

0.9396 

8.0 

0.9993 

3.9 

0.9594 

8.1 

0.9997 

4.0 

0.9766 

8.2 

1.0000 

4.1 

0.9912 

8.3 

1.0003 

8.4 

8.5 
8.6 

8.7 
8.8 

8.9 
9.0 

9-1 
9.2 

9.3 
.  9.4 
.  9.5 

9.6 

9.7 

9.8 

9.9 

10.0 

10.1 

10.2 

10.3 
10.4 

10.5 
10.6 

10.7 
10.8 

10.9 
11.0 
11.1 
11.2 

11.3 
11.4 

11.5 
11.6 
11.7 
11.8 

11.9 
12.0 
12.1 
12.2 

12.3 

12.4 

12.5 


y(t) 


1.0005 
1.0007 
1.0008 
1.0009 

1.0009 
1.0010 
1.0009 
1.0009 
1.0008 
1.0008 

1.0007 
1.0006 

1.0005 

1.0004 

1.0003 

1.0002 
1.0002 
1.0001 
1.0000 
1.0000 
0.9999 

0.9999 
0.9999 

0.9999 

0.9998 
0.9998 
0.9998 
0.9998 

0.9999 
0.9999 
0.9999 
0.9999 
0.9999 
0.9999 
0.9999 
1.0000 
1.0000 
1.0000 
1.0000 
1.0000 
1.0000 
1.0000 
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Thus 

^  =  (1)(1  +  z^^  )(, 3 )  (53) 

^       ^s^^     2  ^   s^+2.6s3+3.4s2+2.6s+1 

The  block  diagram  is  shown  in  Fig.  7(d).     Following  the  same 
procedures  as  in  the  previous  cases  yields 

Zy  [(1+1  .3T)+z(J+-2.6t+3.^t2+1  .3T3-tW)+z2(6-6.8T2-UW)+ 
z3(-it+2.6T+3.^T2-l  .3T3-rtJ)+z^(1-1 .3T)  j 

_  T^z(z3+1lz2+11z+1)(1+z33)  1..  ^^^^ 

^8(1 -z) 

where    W=T  /6 

Therefore  the  recxirrence  relation  is 

^  (Ay„.,^By^,,-K^,.^^Dy,^)     ^  ^E)  (55) 

■^n  p  F    ^ 

x^  =|o,   1,   12,   23,  2^,   24 24,  25.  36.  47. 

48.  48,    I 

A  =  4+2.6T+3.4T^-1.3T^-W  .       . 

B  =  6.8T^-4W-6 

C  =  4-2.6T-3.4T^+l  .3T^-W 

D  =  1 .3T-1 

E  =  i- 
48 


where 


F  =  1+1.3T 

Equation  (55)  yields 
V     =  0 


n  =  0 


y^   =  S/F  •  n  =  1 

y,  =  (Ay^+12E)/F 


n  =  2 
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(c) 


— tj K 2 
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Ikl 
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^3  = 


n  p 


_  Ay3^+By33-K:y33+Dy3^  +36E 
^35  ■  F 


^36  F 


Results  are  presented  in  Table  12  and  Fig.  6, 


n  =  3 


_  kj^+By^ +23E 

F 
Ayn-1-^^n-2-^n-3-^^nV^^^  ^^n^33 


Ay33-^^3,^3,^30^25E  ^  ^- ^ 

^3^  F 


n  =  35 


a  =  36 


n  ^37 


B,  Second  step  of  compensation 

From  Table  7  and  Table  12  it  is  fo\ind, 

T,  =  2.85  sec. 
do 

"^dl  ~  ^'^^  ^®°*  ■    ' 

Then  '   '      ; .  ' 

Replacing  s  bj  0.67s  in  the  transfer  function  of  equation  (53)  yields 


y  =  (1)(L-L£i)( , 1 )      (56) 

^  ^         2     0.202s^-K).782s3+l  .527s2+1 .7^2s+1        , 

Applying  the  same  process  as  in  the  previous  cases  yields 
Zy  [(0.202+0. 391  T)-i^(0. 871  t3+1.527t2-0.782T-0.808-H//) 

+z2(l  .212-3.045t2-HW)+z2(W-0.871t3+1  . 5271^+0. 782T-0. 808) 
zHo. 202-0. 391  T)] 


Table    12.    y( t  )=  ^""^-^^^ 
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l/l 


+2*6s-^+3.^s^+2.6s+l 


y(0)=:0,    T=0.1    sec. 


■). 
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y(t) 


(t) 


y(t) 


0.0 
0.1 
0.2 
0.3 

0.5 
0.6 

0.7 
0.8 

0.9 

1.0 

1.1 

1.2 

1.3 
l.ii 
1.5 

1.5 

1.7 
1.8 

1-9 
2.0 
2.1 
2.2 

2.3 

2.4 

2.5 
2.6 

2.7 
2.8 

2.9 
3.0 
3.1 

3.2 

3.3 
3.4 

3.5 
3.6 

3.7 

3.8 
3.9 

4.0 

4.1 
4.2 

4.3 

4.4 


0.0000 

4.5 

0.5336 

9.0 

1.0410 

0.0000 

4.6 

0.5487 

9.1 

1.0403 

0.0000 

4.7 

0.5v'^39 

9.2 

1.0392 

0.0001 

4.8 

0.5794 

9.3' 

1.0378 

0.0004 

4.9 

0.5951 

9.4 

1.0361 

0.0010 

5.0 

0.6110 

9.5 

1.0342 

0.0019 

5.1 

0.6271 

9.6 

1.0321 

0.0034 

5.2 

0.6435 

9.7 

1.0298 

0.0055 

5.3 

0.6601 

9.8 

1.0274 

0.0084 

5.4 

0.6770 

9.9 

1.0250 

0.0121 

5.5 

0.6^39 

10.0 

1.0225 

0.0167 

5.6 

0.7110 

10.1 

1.0199 

0.0224 

5.7 

0.7282 

10.2 

•  1.0174 

0.0291 

5.8 

0.745^ 

10.3 

1.0149 

0.0369 

5.9 

0.7625 

10.4 

1.0124 

0.0458 

6.0 

0.7796 

10.5 

1.0100 

0.0559 

6.1 

0.7965 

10.6 

1.0077 

0.0671 

6.2 

O.8132 

10.7 

1.0055 

0.0794 

6.3 

0.8296 

10.8 

1.0034 

0.0928 

6.4 

0.8456 

10.9 

1,0015 

0.1071 

6.5 

0.8613 

11.0 

0.9996 

0.1224 

6.6 

0.8765 

11.1 

0.9980 

0.1385 

6.1 

0.8912 

11.2 

0.9965 

0.1553 

6.8 

0.9053 

11.3 

0.9951 

0.1729 

6.9 

O.9I89 

11.4 

0.9939 

0.1909 

7.0 

O.93I8 

11.5 

0.9928 

0.2094 

7.1 

0.9440 

11.6 

0.9920 

0.2283 

7.2 

0.9556 

11.7 

0.9912 

0.2473 

7.3 

0.9664 

11.8 

0.9906 

0.2665 

7.4 

0.9765 

11.9 

0.9902 

0.2856 

7.5 

0.9585 

12.0 

0.9899 

0.3047 

7.6 

0.9943 

12.1 

0.9897 

0.3235 

7.7 

1.0021 

12.2 

0.9896 

0.3419 

7.8 

1.0091 

12.3 

0.9896 

0.3600 

7.9 

1.0154 

12.4 

0.9898 

0.3776 

8.0 

1.0209 

12.5 

0.9900 

0.3947 

8.1 

1.0257 

12.6 

0.9904 

0.4113 

8.2 

1.0298 

12.7 

0.9908 

0.4274 

8.3 

1.0332 

12.8 

0.9912 

0.4431 

8.4 

1.0360 

12.9 

0.9917 

0.4585 

8.5 

1.0382 

13.0 

0.9923 

0.4737 

8.6 

1.0398 

13.1 

0.9929 

0.4887 

8.7 

1.0408 

13.2 

0.9935 

0.5037 

8.8 

1.0413 

13.3 

0.9942 

O.5186 

8.9 

1.04l4 

13.4 

0.9949 
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_  A(z^+11z^+11z+1)(H^^^)  (57) 

The  recurrence  relation  is 

y        Ay.-l-^^3r„,-K:yn-3-^Dyn-/^    wE^  ^^^  ■      (58) 

n  J.  F 

x^  =  [o,  1.  12,  23,  2i^,   ......  2^^.  25.  36.  47,  ^,  48,   ... 

^,  .......   !• 

where  A  =  0.808-K).782T-l  .527T^-0.871t3-W 

B  =  3.054T^-4W-1.212 

C  =  0.808-0. 782T-1.527T^-H). 871  T^-W 

D  =  0.391 T-0. 202 

E  =  T^/48 

W  =  T^/6 

F  =  0.391 T+0. 202  .         . 

Equation  (58)  yields 


•'O 

E 

yi  =F 

Ayi+12E 

'^  F 


Ayg+By-j  +23E 
^3  = i 


Ay22+^721  +Cy20:*M9:*^^ 
^23  = i^ 


n  =  0 

n  =  1 

n  =  2 

n  =  3 


n  F 


n  =  23 


Ayp,-^gy,,^/2l-^Qy20^36E       .•  ^  ^  ^^ 

^2if  V 


4 


^2 

F 
Results  are  presented  in  Table  13  and  Fig.  6.  The  maximum  overshoot 
reduces  from  0.1118  to  0.0401  and  the  time  delay  reduces  from  2.85  seconds 
to  2.83  seconds. 


^J 


Table  13.     y(t) 


=I-v 


(1    +222)/2     (1/s) 


0.202s^+0.782s3+1 .527s2+1 .742s+1 


-).y(o)=o, 


T=0.1    sec, 


y(t) 


y(t) 


y(t) 


0.0 

0.0000 

4.4 

0.8740 

8.8 

0.9920 

0.1 

0.0000 

4.5 

0.8958 

8.9 

0.9930 

0.2 

0.0001 

4.6 

0.9164 

9.0 

0.99^ 

0.3 

0.0006 

4.7 

0.9355 

9.1 

,.      0.9950 

OA 

0.0019 

4.8     . 

0.9531 

9.2 

0.9960 

0.5 

0.0043 

4.9 

0.9691 

9.3 

0.9969 

0.6 

0.0082 

5.0 

0.9834 

9.4 

0.9978 

0.7 

0.0140 

5.1 

0.9960 

9.5 

0.9986 

0.8 

0.0219 

5.2 

1 .0069 

9.6 

0.9994 

0.9 

0.0323 

5.3 

1.0161 

9.7 

1 .0001 

1.0 

0.0451 

5A 

1 .0237 

9.8 

1 .0006 

1.1 

0.0604 

5.5 

1 .0297 

9.9 

1.0011 

1.2 

0.0782 

5.6 

1 .0343                 1 

0.0 

1.0015 

1.3 

0.0984 

5,7 

1 .037^ 

0.1 

1.0018 

\A 

0.1207 

5.8 

1 .0393 

0.2 

1 .0020 

1.5 

0.1450 

5.9 

1 .0401 

0.3 

1.0021 

1.6 

0.1708 

6.0 

1.0398 

0.4 

1 .0022 

1.7 

0.1978 

6.1 

1 .0387 

10.5 

1 .0022 

1.8 

0.2258 

6.2 

1 .0368 

10.6 

1 .0021 

1.9 

0.2542 

6.3 

1 .0342 

10.7 

1.0020 

2.0 

0.2828 

6.4 

1.0312 

10.8 

1.0018 

2.1 

0.3112 

6,5 

1 .0279            •     ■ 

10.9 

1.0016 

2.2 

0.3390 

6.6 

1 .0243 

11.0 

1.0013 

2.3 

0.3658 

6.7 

1 .0205 

11.1 

1.0011 

ZA 

0.3916 

6.8 

1.01 67 

11.2 

1 .0008 

2.5 

0.4164 

6.9 

1.0130 

11.3 

1 .0005 

2.6 

0.i^2 

7.0 

1 .0094 

11.4 

1 .0002 

2.7 

0.4633 

7.1 

1.0059 

11.5 

1 .0000 

2.8 

0.4860 

7.2 

1.0027 

11.6 

0.9997 

2.9 

0.5084 

7.3 

0.9998 

11.7 

0.9995 

3.0 

0.5310 

7.4 

0.9972 

11.8 

0.9992 

3.1 

0.5537 

7.5 

0.9949 

11.9 

0.9990 

3.2 

0.5769 

7.6 

0.9930 

12.0 

0.9989 

3.3 

0.6006 

7.7 

0.9914 

12.1 

0.9987 

3A 

0.6248 

7.8 

0.9902        '     ^■ 

12.2 

0.9986 

3.5 

0.6495 

7.9 

0.9893 

12.3 

0.9985 

3.6 

0.6746 

8.0 

0.9887 

12.4 

.    0.9984 

3.7 

0.7001 

8.1      ■ 

0.9884 

12.5 

0.9984 

3.8 

0.7259 

8.2 

0.9883 

12.6 

0.9983 

3.9 

0.7516 

8.3 

0.9885 

12.7 

0.9983 

4.0 

0.7772 

8.4 

0.9889 

12.8        ■■ 

0.9984 

4.1 

0.8024 

8.5 

0.9895 

12.9 

0.9984 

4.2 

0,8271 

8.6 

0.9902 

13.0 

0.9985 

4.3 

0.8510 

8.7 

0.9911 

13.1 

0.9985 

i|4 
DISCUSSION 

Investigation  of  the  graphs  and  numerical  results  in  the  uncompensated 
cases  shows  that  as  the  order  of  the  Butterworth  filter  increases  the  ratio 
of  rise  time  to  half-cycle  overshoot  duration,  mT,  increases;  this  behavior 
is  shown  in  RLg.  8.  Thus  digital  filter  compensation  is  most  effective  on 
higher  order  Butterworth  filters  and  a  smoother  transient  response  curve 
results  as  seen  in  Flg.'s  4,  5.  6. 

The  transient  response  of  the  higher  order  Butterworth  filter  with 
unit  step  input  always  has  larger  maximum  overshoot  and  time  delay  which 
cannot  be  extensively  reduced,  the  first  time  when  the  digital  filter 
compensation  is  introduced.  However,  the  process  can  be  repeated  until 
the  desired  reduction  in  transient  overshoot  is  obtained;  this  is  the  goal 
of  digital  filter  compensation. 

For  comparison,  posicast  control  compensation  is  applied  to  each  case 
and  the  resulting  transient  responses  are  presented  in  Rig.'s  ^J-,  5.  6.  In 
the  second  order  case,  the  two  compensating  methods  give  almost  the  same 
result,  although  there  is  inflection  in  the  digital  filter  compensation 
case.  In  the  third  order  case  large  undershoot  is  produced  with  posicast 
compensation  and  digital  filter  compensation  gives  better  results.  In  the 
fourth  order  case,  oscillation  with  much  larger  undershoot  occurs  under 
posicast  control  compensation;  however,  digital  filter  compensation 
results  in  much  better  behavior. 
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The  essential  content  of  this  report  is  to  analyze  Butterworth 
filters  in  the  time  domain  by  using  a  digital  method  for  finding  their 
transient  performances.  Digital  filter  compensation  is  effected  by 
multiplication  in  the  Z-transform  domain  by  (l+z^)/2.  For  comparison, 
posicast  compensation  of  a  compensator  (1-P)+Pz^  is  employed  in  each 
digital  filter  compensation  case.  Tna  approximate  digital  analysis  is 
accomplished  by  the  Naumov  program  which  is  an  improved  accuracy  trape- 
zoidal convolution  program. 

Transient  responses  of  higher  order  digital  Butterworth  filters 
employing  digital  filter  compensation  are  much  better  than  those  employing 
posicast  compensation.  Digital  filter  compensation  can  be  repeated  until 
desired  reduction  in  transient  overshoot  is  obtained—the  goal  of  digital 
filter  compensation. 

This  report  introduces  an  interesting  subject,  the  pseudo-Butterworth 

function.  Pseudo-Butterworth  functions  can  be  converted  into  approximate 

l/n 
Buttenxorth  functions  by  replacing  s  with  s/(b^)  '  ,  where  b^  is  the 

coefficient  of  the  n"*^  term  of  pseudo-Butterworth  polynomial  of  order  n. 


